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We study the ground-state density profile of a spin-orbit coupled / = 2 spinor condensate in a 
quasi-one-dimensional trap. The Hamiltonian of the system is invariant under time reversal but not 
under parity. We identify different parity- and time-reversal-symmetry-breaking states. The time- 
reversal-symmetry breaking is possible for degenerate states. A phase separation among densities 
of different components is possible in the domain of time-reversal-symmetry breaking. Different 
types of parity- and time-reversal-symmetry-breaking states are predicted analytically and studied 
numerically. We employ numerical and approximate analytic solutions of a mean-field model in this 
investigation to illustrate our findings. 
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I. INTRODUCTION 

Since the first experimental realization of the spinor 
Bose-Einstein condensate (BEC) with a gas of 23 Na 
atoms trapped in an optical trap |Tj, a lot of theoret¬ 
ical and experimental studies have been done on these 
systems HE]- On the theoretical front, mean-field theo¬ 
ries have been developed both for / = 1 j4[ [5] and / = 2 
spinor Bose-Einstein condensates (BECs) [51 [7. In our 
previous work [5], we studied the ground state structure 
of an / = 1 spin-orbit (SO) coupled spinor BEC with 
fixed magnetization in a quasi-one-dimensional (quasi- 
ID) trap [5] within the framework of the mean-field the¬ 
ory. In the present work, we extend this study to study an 
SO coupled / = 2 spinor condensate in a quasi-ID trap. 
The SO coupling, which relies on the generation of the 
non-Abelian gauge potentials coupling the neutral atoms 
m, can be experimentally realized by controlling the 
atom-light interaction. A variety of SO couplings can be 
engineered by Raman dressing the hyper-fine states. The 
parameters of atom-light interaction Hamiltonian, and 
hence those of coupling can be controlled independently 
m- The SO interaction with equal strengths of Rashba 
[12] and Dresselhaus [HQ3 couplings has been achieved 
recently p5UT5j . The experimentalists employed a pair of 
Raman lasers to create a momentum-sensitive coupling 
between two internal atomic states of 87 Rb [15j [lBj . This 
has lead to a flurry of other experiments on SO-coupled 
pseudo-spinor BECs PU- The generation of SO cou¬ 
pling involving the three hyper-fine spin components of 
an / = 1 spinor condensate using Raman dressing has 
also been studied theoretically m in]- Recently, SO 
coupling has also been experimentally realized in degen¬ 
erate Fermi gases of 40 K and 6 Li [2UI . 

Wang et al. m studied theoretically the ground states 
of pseudo-spin-1/2 two-component BEC with SO cou- 
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pling and of three-component / = 1 spinor BEC. In the 
presence of SO coupling, the ground states of all spinor 
BECs — pseudo-spin-1/2 [22], / = 1 [23], and 2 [23] may 
exhibit different types of nontrivial density distribution. 
In the presence of a uniform magnetic field the ground 
states of / = 1 HH] [25] [25] and / = 2 [2] [7] [27] spinor 
BECs exhibit interesting behavior including the possibil¬ 
ity of a phase separation mi m\ among different spin 
components. There have also been different static and 
dynamic studies in SO-coupled BECs, such as, Joseph- 
son oscillation [25], intrinsic spin-Hall effect [2S], solitons 
130] , force on a moving impurity [31] , chiral confinement 
[32], and superfluidity [33], etc. 

In this paper, we investigate the ground state of an 
SO-coupled / = 2 spinor BEC in a quasi-ID trap for an 
arbitrary magnetization. The Hamiltonian of this system 
is invariant under time reversal T but not under parity. 
Consequently, different types of parity-breaking states 
are found. Time-reversal symmetry-breaking states are 
also found in the presence of degeneracy. In the ab¬ 
sence of degeneracy, the states preserve time-reversal 
symmetry. The five spin-component wave functions of 
the / = 2 spinor condensate satisfy a coupled mean- 
field Gross-Pitaevskii (GP) equation with three interac¬ 
tion parameters: Cq cx ( 4a2 + 3a4)/7, Ci oc (<24 — a2)/7, 
and ci oc (7ao — 10a2 + 3ei4)/7, where a 0, 02, and <24 
are the s-wave scattering lengths in total spin f to t = 0, 2 
and 4 channels. The whole Ci versus C2 parameter space 
can be divided into sub spaces with distinct symmetry 
properties of the densities of spinor components. We 
have the ferromagnetic phase for Ci < 0 and C2 > 20ci, 
anti-ferromagnetic phase for C2 < 0 and C2 < 20ci, and 
cyclic phase for ci > 0 and C2 > 0. In the ferromagnetic 
phase increasing magnetization lowers energy, whereas 
in the anti-ferromagnetic phase the lowest energy is at¬ 
tained for zero magnetization. Miscible configuration for 
five component densities of an SO-coupled spinor BEC is 
obtained for C2 < 20ci and a phase separation is possi¬ 
ble for C2 > 20ci. We find that for sufficiently strong 
spin-orbit coupling, the SO-coupled spinor condensate 
has atoms only in to/ = 2 and nif = —2 states. Time- 
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reversal-symmetry is preserved for states only in the anti¬ 
ferromagnetic phase, and can be broken in other phases. 
As parity is not a good quantum number, it is broken in 
all domains. We use the numerical solution of the gener¬ 
alized mean-held GP equation [33] for this investigation. 


The paper is organized as follows. In Sec. |TTJ we de¬ 
scribe the coupled GP equations used to study the SO- 
coupled / = 2 spinor BEC in a quasi-ID trap. From an 
analytic consideration of energy minimization, we pre¬ 
dict the expected density profiles for different set of pa¬ 
rameters. Specifically, we predict the parameter space, 
where a phase separation can take place. We also predict 


different types of parity-breaking states. In Sec. Ill we 
numerically study the SO-coupled spinor BEC in a quasi- 
ID trap. We identify the different types of symmetry¬ 
breaking states and phase-separated states in different 
parameter domains. We conclude by providing a sum¬ 
mary of this study in Sec. [TV] 


II. MEAN-FIELD MODEL FOR A SO-COUPLED 

BEC 


Spin-orbit coupling can be generated for the hyper- 
fine states of neutral atoms by suitably controlling the 
atom-light interaction. The idea was realized experimen¬ 
tally by Lin et al. |15] for two hyper-spin components of 
the 87 Rb hyper-hne state 5S!/ 2 employing two counter- 
propagating Raman lasers of wavelength A r oriented at 
an angle /3 r . This lead to the SO coupling with strength 
7 = hkr/m, where k r = 27rsin(/3 r /2)/A r and m is the 
mass of an atom. However, here we consider the SO cou¬ 
pling among the five spin components of the / = 2 state, 
e.g., |/ = 2 ,to/ = 2), |/ = 2,m/ = 1), |/ = 2 ,m f = 0), 
|/ = 2,m/ = —1) and |/ = 2, to/ = —2), where to/ is the 
2 projection of /. By generalizing the method discussed 
in Ref. [15], this SO coupling among the five hyper-fine 
spin components can be generated by engineering a suit¬ 
able atom-light interaction Hamiltonian as in Ref. [13- 

In order to realize a quasi-ID SO-coupled spin-2 BEC 
along x axis, we consider a trapping potential with angu¬ 
lar frequencies along y and z axes much larger than that 
along x axis. The resultant strong transverse confine¬ 
ment ensures that the dynamics is frozen along y and z 
axes. Then, the single-particle quasi-ID Hamiltonian of 
the system under the action of a strong transverse trap 
of angular frequencies ui y and u z along y and z axes, 
respectively, can be written as mm 

Ho = + 7 Px^z + ( 1 ) 

2 TO 

where p x = — ihd x is the momentum operator along x 
axis, H is the Rabi frequency mm, and E- and E x are 
the irreducible representations of the z and x components 


of the spin-2 matrix, and are given by 


E 2 


/2 0 0 0 0 \ 

0 10 0 0 

0 0 0 0 0 

000-1 0 

V 0 0 0 0 -2 / 


£ 


X 


(0 1 0 0 0\ 

1 0 y/| 0 0 

0 >Jl 0 0 

0 0 yjl 0 1 

V o o o 10 / 


( 2 ) 


(3) 


A. Mean-field model in a quasi-ID trap 


If the interactions among the atoms in the BEC are 
taken into account, in the Hartree approximation, using 
the single particle model Hamiltonian (JT|) , a quasi-ID [0] 
spin-2 BEC can be described by the following set of five 
coupled mean-field partial differential GP equations for 
the wave-function components ipj [21 G33 




H 2 d 2 


ih ~fo= ( + V{x) + c 0 p) ip 2 ~ 2 ihj— 


2 to dx 2 


dip 2 


dx 

H^i, (4) 


+ ci (F-ipi + 2 F z ip 2 ) + (c 2 /v / 5)0V’* 

dp,i ( h 2 d 2 \ dip i 

= v(x>+Cop ) * - ,hl 
+ c i(\/372-F-V’o + F + ip 2 + F z p >i) - (c 2 /a/5)0^* : 

+ H [i />2 + (v / 3/2) , 0o]; (5) 

dp’ 0 


ih 


dt 


h 2 d 2 
2to dx 2 


C 2 


+ V(x) + c 0 p ) ipo + -^=9$ 


+ ^c 1 (F_V’-i +F+ipi) + + pJ-\), (6) 

dip -1 
dx 


..dip -1 f H 2 d 2 \ 

+ c 1 (\/ r 3j2F + ip 0 + F-ip -2 - F z ip-i) - (c 2 /V$)Qipl 

+ n[(V3/2)ip 0 + ip- 2 ], (7) 


ih 


dp ,-2 ( h 2 d 2 \ dp ,-2 

w = (7 ss + CM+ cop) 

+ C 1 (-F+V’-1 ^ 2 F z ip_ 2 ) + (c 2 /V5)0^2 + ttip-i, (8) 


where V(x) = mui 2 x 2 /2 is the ID harmonic trap, cq = 
2fi 2 (4a 2 + 30 , 4 ) / (7ml 2 z ), Ci = 2 h 2 (a±-a 2 )/{ 7 ml 2 z ), c 2 = 
2h 2 (7ao — 10a 2 + <$( 14 )/(7ml 2 z ), a 0 , a 2 and 04 are the s- 
wave scattering lengths in the total spin / tot = 0, 2 and 4 
channels, respectively, pj = ipj \ 2 with j = 2, 1 , 0, — 1 , —2 
are the component densities, p(x) = J2j =-2 Pi the to¬ 
tal density, and l yz = ^/hj(moj yz ) with uj yz = ^uj y ui z is 













3 


the oscillator length in the transverse y — z plane and 

F+=F*= 2(ip%ipi + ip-id- 2 ) + + V’oV’-i)) 

F z = 2(\H 2 -W-2\ 2 ) + \^\ 2 -W-i\ 2 , 

_ 2lp 2 1p-2 - Zi’llt’-l + V’O 

V5 

Here F = (F x ,F y ,F z ) is the spin density vector; F± = 
F x ± F y . The normalization condition is J_ dxp(x) = 
N , where N is the total number of atoms. In order to 
transform Eqs. Q - Q into dimensionless form, we use 
the scaled variables defined as 

t = ui x t, x=^~, (j)j(x,t) = —^ipj(x,t), (9) 

1 0 v N 


where Iq = y/h/(mu) x ) is the oscillator length along x 
axis. Using these dimensionless variables, the coupled 
mean-field Eqs. ([4]) - ([ 8 j in dimensionless form are 


.d(t> 2 


dt 


1 d 2 


+b + cop U 2 - 217—r 


2 dx 2 


.d<t>2 


dx 


+ ci(F-di + 2 hfo) + (c 2 /V5)e4>*_ 2 + 001 , ( 10 ) 


.d<l> 1 


dt 


1 d 2 


+ V + 5 0 p <f>i — *7^v- 


2 dx 2 


dcj) 1 


dx 


+ ci( yfzj2F-cf>o + F+fo + Fs<t> 1 ) - (d 2 /V 5)0f_i 
+ ft [<f>2 + (\/ 3 / 2 )<^o]> 


.d<t> 0 


dt 


1 d 2 


+^ + coP 


2 dx 2 




+ (F_<j>~i + F + (j> n) + + d-i ), 


( 11 ) 


( 12 ) 


■ 50_i 
1 dt 


1 5 2 

2 dx 2 


+ V + c 0 p) <j>-i +*7 


<95; 


+ ci(y/3/2F + <t>o + F_<j )_2 — F z cf>-i) - (C2/V5)e# 

+ H [ (x/3/2) 0o + <^- 2 ] > (13) 

.d<j>-2 
dt 


1 <9 2 


= ' ^2aj2 +^ + C0P ) ^-2+2*7 




<95; 


+ Ci(F+^_i - 2F5</>_ 2 ) + (c 2 /v / 5)0</>2 +U^_i, (14) 


where V = i 2 / 2 , 7 = hk r /(muj x lo), fl = fl/(hui x ), 
c 0 = 27V(4a 2 + 3a 4 )/(7Z 2 J, cj. = 27V(a 4 - a 2 )/(7Z 2 J, 
c 2 = 2N(7a 0 - 10a 2 + 3a 4 )/(7Z 2 J, pj = \(f)j\ 2 with 
J = 2 , 1 , 0 ,- 1 ,- 2 , and p = Ej =- 2 l<A ?| 2 and 


F+ = F* =2(cf>* 2 di + £.^- 2 ) + \/ 6 (^*^o + ^50_i), 

Fz = 2(|</> 2 | 2 — | 0 — 2 1 2 ) + | 0 i | 2 — | 0 —i | 2 , 

- _2<t>2<t>-2 - 20 1 0_ 1 + </>§ 

%/5 

The normalization condition satisfied by dj 's is 
j^ > oo p(x)dx = 1. One of the aims in the present work 


is to find the ground state of an / = 2 spinor condensate 
with a fixed magnetization, which is defined by 

/ OO 

F z dx. (15) 

-OO 

Depending on the values of c 4 and c 2 the system in the 
absence of magnetic field and SO coupling can have a 
variety of ground states [2j. For the sake of simplicity of 
notations, we will represent the dimensionless variables 
without tilde in the rest of the paper. 


B. Uniform BEC: Analytic Consideration 


The energy of a uniform (trapless) spinor BEC in the 
presence of SO coupling and magnetic field is [21 SU 


E = N 


{ I 

i =-2 


d(j)j 


dx 


J = - 2 


+ (co P~ + ci | F | 2 + c 2 | 0| 2 )/2 + ft [<)>2 di + di 02 
+ d-ld-2 + $L 2 </>-l + ('v/3/2) (dido + < / ) 0 < / > l + 
00 ^-1+^-1^0)] ( 16 ) 


First we consider the formation of spatially-separated 
non-overlapping (phase-separated) states from a consid¬ 
eration of energy minimization. As in the case of a 
/ = 1 spinor BEC 0, the energy term proportional to 
cq in Eq. (16) can not lead to a phase separation as 


it contains terms Nc 0 J(p 2 /2 + p' 2 /2 + pjpt)dx , where 
j,j' = 2 , 1 , 0 , 1 , —2 and j ^ /, and hence corresponds to 
a situation where inter- and intra-species interactions are 
of equal strength. The situation is analogous to a binary 
BEC with a ai 2 = y / ana 22 , where an and a 22 are intra¬ 
species and ai 2 the inter-species scattering lengths. Such 
a binary BEC has equal strengths of inter- and intra¬ 
species nonlinearities and is always miscible in the pres¬ 
ence of a ID harmonic trap [361137]. The interaction 
energy of the / = 2 spinor condensate in the absence 
of SO coupling and magnetic field (7 = fl = 0) can be 
written as 


Cl 


E;„ t =N M -i 


2 L 


4p 2 + Pi + P-i + 4 p 2 _2 

+ f>p 0 (pi + p- 1 ) + %y/pip2P-\P-2 cos(6Li_ 2 - 0 2i ) 

+ 12poV P1P-1 cos(0 O -i — 0io) + 4 v / 6 ^p-i- v /p- 2 Po 
X cos(0_i _ 2 - 0 O -l) + y/p-ip-2p0pl cos(0_i _ 2 - 0io) 

+ y/P2P1P0P-I cos(0 2 i — 0 o_i) + Ply/P2P0 cos(0 2 i — 0io)) 

+ 8p 2 pi + 8p_ip_ 2 - 4 p 2 p_i - 8p 2 p_ 2 - 2pip_i 


- dpip—2 


+ ^ 

10 L 


4 p 2 p—2 + 4 pip_i - 8y/p 2 piP-2P-i 
X cos( 0 2 i — 0_1_ 2 ) + 4 p 0 (y/P2P-2 COs( 0 2 o — 00-2) 
-\/PlP-l COS( 01 Q - 0 O _i)) + Pq } da:: > (l 7 ) 
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where the wave-function component <pj is written as (pj = 
^/Jjj exjp(i 6 j) with 9j the phase. The phase difference 
between the zth and j th components is written as 9ij = 
9i — 9j. 

To understand the phase separation and spontaneous 
symmetry breaking of the states, we consider the sta¬ 
tionary eigenvalue problem of the lowest-energy state of 
a uniform non-interacting system with SO coupling while 
Eqs. (fTol) - (fl4l) become 


Ecpj ( x ) = N 


1 8 2 
2 dx 2 


*77 


d_ 

dx 


(f>j{x)- ( 18 ) 


The two independent solutions of Eq. (18) for the lowest- 
energy state are <p ±2 = a ±2 exp(=F 2 z 7 £Jwith normaliza¬ 
tion |a 2 | 2 + |cv_ 21 2 = 1 and <j)± i = 0 ± o = 0 with en¬ 
ergy E = E m = —2N"/ 2 . The components j = ±1,0 
have higher energies. In the presence of a trap and inter¬ 
actions, in the actual numerical calculation, for a suffi¬ 
ciently large SO coupling 7 only the components j = ±2 
survive. Hence, the analytic solutions of Eq. (18) is 


very useful to understand many features of the actual 
numerical solution. It is clear from Eq. |l 8 | that these 
plane-wave solutions will lead to smooth density profiles 
in the presence of a trap while their real and imaginary 
parts will, in general, have oscillating behavior. In the 
case, when the two components j = ±2 survive, the in¬ 
teraction energy becomes 


£int = 2 TV 


Ci(p 2 + P-2) 2 + 


c 2 — 20 ci 


-P 2 P -2 


dx. 

(19) 


In Eq. (191, only the product term P 2 P -2 controls the 
phase separation between components ±2. A repulsive 
(positive) product term will facilitate a phase separation, 
as the energy can then be minimized by reducing the 
overlap between the components. This will happen for 
c 2 > 20 ci. 


In the presence of SO coupling the Hamiltonian is 
invariant under time reversal T but not under parity. 
Hence, parity is not a good quantum number. Different 
types of simple parity-breaking states are found. The 
non-degenerate states should possess time-reversal sym¬ 
metry. However, a pair of degenerate states, which trans¬ 
form into each other when operated upon by T, break 
time-reversal symmetry. For a sufficiently large 7 , when 
only the components j = ±2 survive, the time rever¬ 
sal symmetry is broken for the phase-separated profiles, 
whereas the miscible profiles preserve time-reversal sym¬ 
metry. Next we consider different types of symmetry¬ 
breaking states. 

First we consider different overlapping parity-breaking, 
nevertheless time-reversal-symmetric, states. The parity 
property of <f>± 2 and the relations between the real and 
imaginary parts, denoted by 7Z and X respectively, of 
<p 2 and (f >-2 for some of these parity-breaking and time- 
reversal-symmetric states are listed in Table [TJ In these 
examples the real and imaginary parts of the wave func¬ 
tion may have definite parity, but not the total wave 


function, as parity is not a good quantum number. How¬ 
ever, no such simple relation is obtained for a general 
a± 2 , where neither the real part nor the imaginary part 
of the wave-function components have a definite parity. 
Also, similar symmetry-breaking states are expected for 
the j = ±1 component states when they are nonzero. 
These types of symmetry-breaking states were found in 
the actual numerical calculation in the presence of trap 
and interaction terms. 

Now we consider some examples of time-reversal 
symmetry-breaking states in the presence of SO coupling. 
These states are phase-separated (non-overlapping). 
There could be a complete phase separation between the 
j = ±2 components when the two components symmet¬ 
rically move to two sides of x = 0. In that case, sup¬ 
pose the components j = ±2 are centered at x = ±Xo, 
then there will be no definite parity of the real and 
imaginary parts, but one can have properties, such as, 
TZ[(j) 2 {x - £ 0 )] = ±TZ[(f>- 2 {x + a; 0 )], X[0 2 (z - Zo)] = 
±X [4>- 2 {x + zo)], or TZ[(t> 2 {x - x 0 )\ = ±X[0_ 2 (® + x 0 )], 
TZ[(p- 2 {x + £ 0 )] = ±X[ 0 2 (^ — £ 0 )], etc- In these cases the 
densities break the symmetry of the trapping potential: 
Pj(x) ^ pj(-x). There could be another type of phase 
separation where one of the components, say j = — 2 , 
stays at the middle and the other component j = 2 
breaks into two parts and stay symmetrically on both 
sides of origin. In this case, the real and imaginary parts 
of the middle component (j = — 2 ) have opposite parities, 
and the real and imaginary parts of the outer component 
(j = 2 ) either map into each other or have opposite pari¬ 
ties: 7 Z[ 4 > 2 {x)\ = ±l[(j) 2 {—x)\, or lZ[<f> 2 (x)\ = ±lZ[<f> 2 (—x)] 
and X[<^ 2 (x)] = q=X[</> 2 (—x)]. In these cases, the symme¬ 
try of the trapping potential is reflected in the densities: 

pA x ) = Pj ( x ) ■ 


For a moderate SO coupling, in the trapped system 
<)>o = 0 and interesting conclusions can be reached ana¬ 
lytically in such a case. In the limit |0o| 0, the inter¬ 

action energy is 


Ei n t — N J | — [4(p 2 + pi + p-i + p- 2) 2 — 3(pi + p- 1) 2 

+ %\JP 2 P- 2 P\P-\ cos(0_i_2 — 0 2 1 ) — 12p 2 p~l 
- 16p 2 P-2 - 4pip_i - 12pip_ 2 ] + ^ [4p 2/ o_2 

± 4 P 1 P -1 — P 2 P- 2 P 1 P -1 cos(6» 2 i — 0_i_ 2 )] jcte. (20) 


The positive terms in Eq. (20 1 involving a product of 
different density components should enhance a phase sep¬ 
aration, as they can be minimized by making the overlap 
between the component wave functions zero. For the ex¬ 
tremum values of cos(0 _!_ 2 — 9 2 i) = ±1, the energy of 
Eq. (20 1 can be written as 
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O' 2 

OL — 2 

Parity property of <j >±2 

Relation between <j>2 and 0 2 

±1 

x/2 

±1 

y /2 

7 £[ 0 ±2(2:)] = U[<j)±2(-x)], I[0±2 (*)] = -X[<t>±2(-x)] 

ll[(j>2(x)] = 7 ^[ 0 _ 2 (a:)], X\<j>2{x)} = -X[<j>-2(x)\ 

±1 

\/2 

V 2 

Tl[<t>±2{x)] = 7£[0±2(-£)], X[0±2(2)] = -X[<j>±2(-x)] 

TZ[ 4 > 2 (x)] = -H[<f)-2(x)\, X[(/)2(x)\ =X[<j)-2(x)] 

V 2 

±£ 

V 2 

2(x)\ = -TZ[(/>±2 {-x)], I[(f>± 2 {x)\ = X[<j>±2(-x)\ 

ll[(j)2(x)] = -U[<f>-2(x)}, X[<j)2(x)] =I[<j)-2{x)\ 

\/2 

ii 

y /2 

1 l[<t>± 2 {x)\ = -U[(j>±2{-x)\, X[<j>±2(x)\ = X[ 4 >± 2 {-x)\ 

TZ[ 4 > 2 (x)] = 71 [ 4 >- 2 (x)\, X\<j>2{x)} = -X[(j>-2(x)\ 

i 

V 2 

±1 

\/2 

H[(p2{x)\ = -U[(t)2 (-*)], TZ[c/)-2(x)] = 7 Z[(j>-2(-x)], 
X[<t>2(x)\ =X[(j> 2 {-x)], X\(j>-2{x)\ = -X[tj>- 2 (-i r)j 

TZ[(/)2(x)] = ±X[(/)-2(x)}, T[ 0 2 (*)] = ± 7 Z[(j>-2{x)) 

±1 

V 2 

i 

7! 

lZ[(j)2{x)] = 1 Z[ 4 > 2 (—x)\, TZ[(j)-2{x)\ = -TZ[(j>-2 (-*)], 
X[<t> 2(2:)] = -X[<j> 2 (-x)], X[<f>- 2 (x)\ = X[tf>- 2 (-x)] 

7 Z[<j>2{x)\ = ±X[ 4 >- 2 {x)\, X[(/)2{x)\ = ± 1 Z[(f)- 2 {x)\ 


TABLE I: Parity-breaking states with different choices of a ±2 . The third column defines the parity of the real and imaginary 
parts of 0± 2 , whereas the fourth column shows the relations between the real and imaginary parts of <f >2 and <f>- 2 . 


Ei n t — — J |4ci(p2 + Pi + P-i + P- 2) 2 ~ 12cip2P-i 

4 _ _ 2 

— 12 ci/?i/j _2 + g( c 2 20ci) P 2 P -2 T \/PiP-i] 


- 3ci(pi - p- 1) 2 T 2Aci^p 2 p-2PiP-i}dx. 


( 21 ) 


Let us in addition consider zero magnetization: f dxpj = 
f dxp-j,j ^ 0. For ci > 0, C 2 > 20ci, the crossed terms 
in density involving pj and p-j are positive and those 
involving p± 2 and p^\ are negative. Hence a stable state 
with energy minimization will correspond to a phase sep¬ 
aration between components ±2 and between ±1 while 
maintaining overlap between components 2 and —1 and 
between —2 and 1. In case of this phase separation the 
term ^2Ac\J~p2pZ2p\p^\ will contribute zero. Similarly, 
for ci > 0 and C 2 < 20ci, the dominating contribution of 
the terms in Eq. (211 involving a product of densities are 


negative, and can be minimized by increasing the over¬ 
lap between components and one can never have a phase 
separation. For Ci < 0 and C 2 > 20ci, the dominating 
contribution of the terms in Eq. (211 involving a prod¬ 


uct of densities are positive, and can be minimized by 
accommodating all the atoms in phase-separated m/ = 2 
and —2 components. All the atoms in m/ = 2 and 
—2 components also ensure the minimum contribution 
from the — 3Nci f(pi — p^i) 2 /2dx term. Finally, for 
Ci < 0, C 2 < 20ci, the crossed terms in density involving 
Pj and p-j are negative and those involving p ±2 and p^i 
are positive. Hence if a phase separation occurs it will 
be between components 2 and —1 and between —2 and 1 
while maintaining overlap between components ±2 and 
between ±1. However, a consideration of minimization 
of the repulsive contribution —3Nc\ J (pi — p-\) 2 /2dx to 
energy (211 for overlapping ±pj requires pi = p_ 1 = 0, 
which will exclude the possibility of a phase separation. 

There are several known phases of this spin-2 sys¬ 
tem. For ci < 0, C 2 > 0, the state of largest magneti¬ 
zation corresponds to the lowest-energy state and such 
states are termed ferromagnetic. Even in the absence 


of SO coupling these states violate time-reversal sym¬ 
metry. For Ci > 0, C 2 < 0, the state of zero magneti¬ 
zation has the lowest energy corresponding to the anti¬ 
ferromagnetic, or polar, or nematic phase. These states 
with M = 0 preserve time-reversal symmetry in the ab¬ 
sence as well as presence of SO coupling. The ferromag¬ 
netic and anti-ferromagnetic phases also extend to the 
domain — Ci < 0, C 2 < 0 — however, separated by the 
line C 2 = 20ci. For Ci > 0, C 2 > 0 neither ferromagnetic 
nor anti-ferromagnetic property prevails and a new phase 
termed cyclic emerges. A separation of phase is expected 
for ferromagnetic material and not for anti-ferromagnetic 
material. In this domain of cyclic phase, the time-reversal 
symmetry is broken in the absence of SO coupling. In 
the presence of a sufficiently strong SO coupling, when 
only the components j = ±2 survive, the time-reversal 
symmetry is broken in the phase-separated domain, i.e., 
C 2 > 20ci, while it is preserved in the miscible domain, 
i.e., C 2 < 20ci. It is interesting that the present analytic 
discussion from a consideration of energy minimization 
could predict the line C 2 = 20ci separating the ferro¬ 
magnetic and anti-ferromagnetic phases corresponding to 
phase-separated and overlapping states, respectively. 

In the following section, by numerically solving the 
coupled Eqs. (10) - (14), we will show that for C 2 > 20ci, 

10 o I 0 is sufficient but not necessary condition for a 

phase separation. 


III. NUMERICAL SOLUTION OF THE 
COUPLED GP EQUATION 

We study the ground state structure of the spinor BEC 
by solving the coupled Eqs. (10) - (141 numerically using 


a split-time-step Crank-Nicolson method 34) [38]. The 
spatial and time steps employed in the present work are 
8x = 0.05 and St = 0.000125. In order to find the ground 
state, we employ imaginary-time propagation. The imag¬ 
inary time propagation neither conserves the normaliza¬ 
tion nor the magnetization as the (imaginary) time evo¬ 
lution operator is not unitary. To fix the normalization, 
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and consequently preserve magnetization, we suggest the 
following approach |39lj . 


A. Calculation of the normalization constants 


The minimization of energy given by Eq. © under 
the constraints of fixed normalization (= 1 ) and magneti¬ 
zation (A4) can be implemented by minimizing the func¬ 
tional 


K = E — fi (j | (f)j 1 2 dx — 1^ — A (^J F z dx — , 

( 22 ) 

where fj, and A are the Lagrangian multipliers and are 
functions o f 0, ’s. The imaginary time equivalent of 
Eqs. (10 1 - (14) using this functional can be written as 


94>j{x,r) 

dr 


SE 




- (H +j\)<j)Ax,T). (23) 


This coupled set of equations is also termed as continuous 
normalized gradient flow equations ||4~0] for / = 2 spinor 
BEC. Applying the first order time splitting to Eqs. (23), 
we break up this equation into two parts 


dr 

d<Pj(x,r) 


5E 


S<t>*Ax,r)’ 


= ~{V + jX)<t>j{x,T), 


(24) 


(25) 


dr \r~ • J-'J 3 

which have to be solved one after the other. The solution 
of Eq. (25) at r = r + St is analytically known: 


<t>j(x, r + St) =dj4>j(x,T) 


= exp 


pT+Sr 


(m + j^) dr 


<t>j{x,r), (26) 


Using this definition of dj , one can derive the following 
relations [55] : 


did_i 

= d 2 , 

(27) 

d 2 d- 2 

= d 2 , 

(28) 

d 2 d 2 _ ± 

= d 3 . 

(29) 


Now, the constraints on norm and magnetization can be 
written in terms of the normalization constants Nj of the 
wave-function components as 

d 2 2 N 2 + d\N x + d 2 0 N 0 + d\N_ x + d 2 _ 2 N_ 2 = 1, (30) 

2 d\N 2 + d\N i - - 2 d 2 _ 2 N_ 2 = M, (31) 

where Nj = f \cj)j(x,T)\ 2 dx. Equations ( [27] ) - ( [3l| ) lead 
to the following set of non-linear algebraic equations 

u 4 N 2 + vu 3 Ni + v 2 u 2 N 0 + v 3 uN_ i + v 4 N_ 2 = 1, (32) 
2u 4 N 2 + vu 3 N i - v 3 uN_i - 2 v 4 N_ 2 = M, (33) 



FIG. 1: (Color online) The c 2 — c\ phase plot illustrating fer¬ 
romagnetic, anti-ferromagnetic, and cyclic phases. The c 2 = 
20 ci line separating the ferromagnetic and anti-ferromagnetic 
phases as obtained from the present analytic consideration is 
shown. Separated phase is possible above this line and misci¬ 
ble phase below this line. 



x// 0 x/Z 0 



FIG. 2: (Color online) Component densities of a 23 Na spinor 
BEC of 10000 atoms for 2 = 0 and Co = 242.97, ci = 12.06 > 
0, and c 2 = —13.03 < 20ci. The parameters used are (a) 
7 = A4 = 0, (b) 7 = M. = 0, (c) 7 = 7 ', M = 0, and (d) 
7 = 7 ', At = 0.5; here 7 ' can have any arbitrary real value 
including 0, and hence denotes any arbitrary strength of SO 
coupling. Time-reversal symmetry is broken in (d). In this 
and following similar figures all quantities are dimensionless 
and only the nonzero density components are shown. 


where u = d\ and v = d%. We use Newton-Raphson 
method m for non-linear system of equations to solve 
Eqs. (32) - (33) after each iteration in imaginary time to 
determine di and do and hence the remaining projection 
operators using Eqs. (27) - (29). 
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B. Numerical Results 


We consider 10000 atoms of 23 Na (in / = 2 hyper-fine 
spin state) in a trapping potential with oj x /(2tt) = 20 
Hz and u) v /(2tt) = uj z /(2tt) = 400 Hz. The oscillator 
lengths with this set of parameters are Iq = 4.69 pm 
and l yz = 1.05 /im. In numerical calculation it is found 
that all over the ferromagnetic domain of Fig. [T] the 
density profiles are qualitatively the same always leading 
to a phase separation. In the anti-ferromagnetic domain 
the density profiles also remain similar without a phase 
separation. In the cyclic domain the density profiles are 
different for C 2 > 20 ci and C 2 < 20 ci, the former leading 
to a phase separation and the latter leading to a miscible 
configuration. We will consider these distinct domains 
in the presentation of results. In addition, in different 
domains we have different types of symmetry-breaking 
states, which we will also illustrate. 

The background s-wave scattering lengths of 23 Na in 
total spin / to t = 0,2, and 4 channels are oo = 34.9 as, 
a.2 = 45.8 ob, and 04 = 64.5as (am, respectively, where 
as is the Bohr radius. With these values of scatter¬ 
ing lengths, we have Co = 242.97, Ci = 12.06 > 0, 
and C 2 = —13.03 < 20ci, corresponding to the anti¬ 
ferromagnetic domain in Fig. |T| In the absence of SO 
and Rabi couplings, and magnetization, 7 = Q = M. = 0, 
there are more than one degenerate ground states [2: . In 
Fig. 2] (a) - (c), three such degenerate ground states 
are illustrated. These states are obtained with differ¬ 
ent initial choices for the wave-function components in 
imaginary-time propagation. For these states, some of 
the wave-function components have zero values. One can 
also have a state where the two components j = ±1 are 
populated with identical density (not illustrated here) as 
in Fig. [ 2 ] (c), which can be generated by a suitable ro¬ 
tation of the state of Fig. [ 2 ] (c) in spin space [2j. In the 
presence of a non-zero SO coupling (7 7 ^ 0, H = M =0), 
the degeneracy between the various ground states is re¬ 
moved and only the non-degenerate state of Fig. [ 2 ] (c) 
with the components j = ±2 survives. In this case there 
is no phase separation as concluded analytically in Sec. 


IIB The degeneracy is also removed for non-zero mag¬ 
netization (7 = 12 = 0,A4 / 0). In this case too, there 
is only one ground state involving components j = ± 2 , 
whose density profile does not change with the introduc¬ 
tion of SO coupling as is shown in Fig. [ 2 ] (d). Hence, 
the degenerate ground states exist only for zero magne¬ 
tization in the absence of SO coupling. 

To obtain Ci > 0 and 0 < C 2 < 20ci, we consider 
ao = 52.35as, 0,2 = 45.8as, and <24 = 64.5ob, leading to 
c 0 = 242.97, ci = 12.06 > 0, and c 2 = 65.76 < 20ci- 
This corresponds to the cyclic domain with miscible 
states. For M = 0, the ground state solution in the 
absence of the SO coupling (7 = 0) is shown in Fig. [ 3 ] 
(a), which can be written as (<^ 2 ,</>i, 4>0i <^- 1 , ( t>- 2 ) T = 
y/p(l, 0, i\/2, 0, 1) t / 2, where T stands for transpose and 
p is the total density. This state has a purely imaginary 
<j)o, a real <f >±2 and <f>±i =0, and is degenerate with the 



FIG. 3: (Color online) Same as in Fig. [ 2 ] for c 0 = 242.97, ci = 
12.06 > 0, and C 2 = 65.76 < 20ci.. The parameters used are 
(a) 7 = M = 0, (b) 7 = 0,Af = 0, (c) 7 = 1,A4 = 0, and 
(d) 7 = 1, At = 0.5. Time-reversal symmetry is broken in all 
cases. 


state = yfp( 1,0,0, v / 2,0) T /v / 3 

shown in Fig. [ 3 ] (b), where the symmetry of the 
states is broken. The latter state has only j = 2,-1 
components. The aforementioned two states break the 
time-reversal symmetry and are degenerate with their 
time-reversed counterparts. For a nonzero 7 , the de¬ 
generacy between these two states is no longer ensured. 
With the introduction of a progressively increasing SO 
coupling 7 , in the state of Fig. [3] (a) p 0 decreases with 
the corresponding increase in p±2 as is shown in Fig. [3] 
(c). With further increase in 7 , po becomes zero and 
only the components j = ±2 survive. The introduction 
of a nonzero magnetization A4 only introduces a split¬ 
ting in the j = ±2 components as shown in Fig. 0 (d)- 
There is no phase separation in this case also. Above a 
critical value of SO coupling 7 , the condensate consists 
of atoms in only to/ = ±2 states, which for A4 = 0 is 
time-reversal symmetric. 

The above study shows that there cannot be a phase 
separation for C 2 < 20ci. Next we consider C 2 > 20ci. 
First we consider a ferromagnetic state with Ci < 0 
and C2 > 0 obtained by employing ao = 52.35ob, 
02 = 45.8as, 04 = 43.0as, leading to cq = 201.36,ci = 
— 1.81, C 2 = 24.15 > 20ci. In this case the densities of the 
spinor BEC of 10000 atoms are shown for 12 = 0 for dif¬ 
ferent values of 7 and A4. in Fig. [4] (a) - (f). From these 
plots we see that with an increase of SO coupling from 
7 = 0 , the overlapping component states separate and 
the population of the j = ±2 components increase at the 
cost of a reduction in population of the j = 0 , ±1 compo¬ 
nents. Finally, for 7 ss> 1 only the components j = ±2 
survive. In all cases a finite non-zero magnetization A4 
breaks the symmetry between densities of the compo¬ 
nents j = ±1 and between j = ±2. The ground state 
solutions in this case are phase separated for 7 greater 









x/Zo 


x/Z 0 



x/Z 0 


x/Z 0 



FIG. 4: (Color online) Same as in Fig. [2] for e 0 = 201.36, ci = 
—1.81,ca = 24.15 > 20ci. The parameters used are (a) 7 = 
M = 0, (b) 7 = 0,M = 0.5, (c) 7 = 0.25, M = 0, (d) 
7 = 0.25, M = 0.5, (e) 7 = 1 , M = 0, and (f) 7 = 1 , M = 0.5. 
Time-reversal symmetry is broken in (b) - (f). 


FIG. 5: (Color online) Same as in Fig. [ 2 ] for Co = 242.97 > 
0,ci = 12.06 > 0 and ci = 459.68 > 20ci. The parameters 
used are (a) 7 = 0.5, M. = 0, (b) 7 = 0.5, M. = 0.5, (c) 
7 = 2,A4 = 0, (d) 7 = 2,A4 = 0.5, (e) 7 = 4, M. = 0, and 
(f) 7 = 4 , M = 0.5. Time-reversal symmetry is broken in all 
cases. 


than a critical value in agreement with the discussion in 
Sec. mm 

Next we consider a cyclic phase with C 2 > 20ci. For 
this we consider ao = 139.6ao, 02 = 45.8ao, and 04 = 
64.5ao leading to cq = 242.97 > 0, Ci = 12.06 > 0 and 
c 2 = 459.68 > 20ci, corresponding to a cyclic phase. As 
the strength of SO coupling 7 is increased, there is a 
phase separation as is shown in Figs. HI (c) -(f)- The 
nature of the phase separation in this case is different 
from that discussed in Ci < 0 and c 2 > 20 ci in the sense 
that the components j = 2 and —1 as well as —2 and 1 are 
overlapping in this case, consistent with the conclusion of 
the theoretical analysis in Sec. |IIB| The solutions in this 
domain always break time-reversal symmetry irrespective 
of the value of 7 . 

Symmetry-preserving versus symmetry-breaking solu¬ 
tions: The imaginary-time propagation, that we use in 
calculation, preserves the symmetry of the initial input 
wave function. Different types of states can be obtained 
with different inputs. For example, the states illustrated 
so far in Figs. [2] — [5]were obtained with Gaussian inputs, 
which were slightly shifted from the origin, for the com¬ 
ponent wave functions. We can obtain different density 
distribution for the components by using initial Gaus¬ 
sian inputs centered at the origin. This is illustrated 
in Fig. [ 6 ] for the same set of parameters as in Fig. [4] 


Co = 201.36, Ci = —1.81, C 2 = 24.15 > 20ci. Distinct 
from Fig. [4j in Fig. [ 6 j the phase separation occurs in 
a different fashion preserving the symmetry of the trap: 
Pi (x) = Pj(-x). One of the components leaves the cen¬ 
tral region and stays symmetrically on both sides of the 
origin. The symmetry-breaking states of Fig. [4] have 
lower energy than the symmetry-preserving states of Fig. 
[ 6 ] with the same sets of parameters. Similar symmetry¬ 
breaking states were found in scalar binary condensates 
|42l. In symmetry-preserving and symmetry-breaking 
cases the phase separation may start with different values 
of magnetization A4 and SO coupling 7 . The magnetiza¬ 
tion and SO couplings in Figs. [ 6 ] (b) and (c) are identical 
to those in Figs. [4] (e) and (d), respectively. For the same 
set of parameters, there is a phase separation in Fig. [4] 
(d) and not in Fig. | 6 ](c). In Fig. [ 6 ] the phase separation 
starts at a larger value of SO coupling. 

Next we show in Fig. [7] some of the possible parity¬ 
breaking states obtained with the parameters of Fig. [2] 
(c), e.g., c 0 = 242.97,ci = 12.06 > 0,c 2 = -13.03 < 
20ci. Figures [7] (a) and (b) were calculated with real 
Gaussian inputs for both the j = ±2 components. 
This will correspond to 0:2 = a _2 = l/x/2 in the dis¬ 
cussion in Sec. EH leading to symmetry properties 
T^[<t>± 2 {x)\ = Tl[<t>± 2 (-a:)], X[<j>± 2 (x)\ = —1[4>±2 (—x)\, 
Tl[<j> 2 (x)\ = TZ[^- 2 (x)], I[<f> 2 {x)\ = —I\4>- 2 {x)\ 1 as il- 
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FIG. 6: (Color online) Same as in Fig. [ 2 ] for c 0 = 201.36, ci = 
—1.81, C2 = 24.15 > 20ci. The parameters used are (a) 7 = 
0.5, M = 0, (b) 7 = 1,M = 0, (c) 7 = 0.25, M = 0.5, and 
(d) 7 = 0.5, A4 = 0.5. Time-reversal symmetry is broken in 
all cases. 




FIG. 7: Real and imaginary parts of the j = ±2 wave- 
function components for the parameters of Fig. [ 2 ] (c). Plots 
(a) and (b) were obtained with two real Gaussian input states 
for j = ±2 components. Plots (c) and (d) were obtained with 
an imaginary Gaussian input for j = 2 and a real Gaussian 
input for j = — 2 components. 


lustrated in Fig. [7] (a) and (b). Figures [7] (c) and (d) 
were calculated with an imaginary Gaussian input for 
j = 2 and a real Gaussian input for j = —2 components. 
This will correspond to a 2 = i/V 2 and a- 2 = l/v^ in 
the discussion in Sec. |II B[ leading to symmetry proper¬ 
ties n[(j) 2 {x)\ = - u[(j) 2 {-x )\, n[(j)- 2 {x)\ = n[(j)- 2 {-x)}, 
Z[fa(x)\ = l[<j> 2 {-x)\, l[<j)- 2 (x)\ = —I[<f>- 2 (—x)], 
TZ[(j) 2 (x)\ = 1[4>- 2 (x)\, X[<j> 2 (x)\ = 1Z[(f>- 2 (x)], as illus¬ 
trated in Fig. [7] (c) and (d). In the same fashion all 
the parity-breaking states obtained in Sec. m from an 
analytic consideration can be realized numerically. 

As in the case of a binary condensate [55] . phase- 


FIG. 8 : Dynamics of a phase-separated spinor condensate 
in the presence of both SO coupling and Rabi term (7 = 
1.0, ft = 0.5 and a zero initial magnetization, (a) The density 
profiles at t = 0, (b) the densities at a i x t = 70.0, and (c) the 
densities at w x t = 275, the time after which the condensate 
has recovered its initial shape after one complete oscillation, 
(d) The component population Nj = J pj(x)dx versus time 
during a complete cycle of the periodic oscillation. 


separated SO-coupled spinor condensates can be catego¬ 
rized either as weakly segregated or strongly segregated. 
In the weakly segregated domain, the total density profile 
preserves the symmetry of the trapping potential (not il¬ 
lustrated here) and has an approximate smooth Gaussian 
profile like that of a single-component BEC in a trap. In 
the strongly segregated domain, a notch appears in the 
total density profile corresponding to symmetry breaking 
solutions shown in Figs. 13(e), (f).| (e) and 11(f)- In 
case of a non-zero magnetization, asymmetrically located 
notch ensures that the total density profiles correspond¬ 
ing to Figs. E](f) and[5](f) do not not have the symmetry 
of the trap (not shown here). 

Spin-mixing dynamics in a phase-separated spinor con¬ 
densate: In the presence of a Rabi term (Q ^ 0) the 
solutions of Eqs. (10 1 - (14), in general, are not station¬ 
ary and exhibit oscillating spin-mixing dynamics |34l I43| . 
To study the spin-mixing dynamics in a phase-separated 
spinor condensate, we again consider 10000 atoms of 23 Na 
with ao = 52.35as, a 2 = 45.8 as a 4 = 43.0as, as in Fig. 
[4j yielding Co = 201.36, ci = —1.81, C 2 = 24.15 > 20ci. 
We first solve Eqs. ( 10 ) - (fl4l) using imaginary-time prop¬ 


agation employing 7 = 1,M. = 0 and f 1 = 0.5 with 
the aforementioned set of parameters. The initial com¬ 
ponent densities so obtained are shown in Fig. [ 8 ] (a). 
We then evolve this solution using real-time propagation. 
The presence of the Rabi term leads to spin-mixing be¬ 
tween the various components. The magnetization M is 
no longer a conserved parameter in the presence of the 
Rabi term as can be interpreted from Fig. H(b), where 
the magnetization is 0.37 at time ui x t = 70. During the 
evolution, the condensate densities execute oscillation pe- 
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riodically recovering the initial shape as shown in Fig. [8] 
(c) at time w x t = 275 after one complete oscillation. In 
Fig. |8](d) we plot the component normalizations Nj ver¬ 
sus time during this periodic oscillation. 

IV. SUMMARY OF RESULTS 

We studied the density profile of a trapped SO-coupled 
/ = 2 BEC for different values of the inter-atomic scatter¬ 
ing lengths. Such modification of the scattering lengths 
can be realized in a laboratory using the Feshbach reso¬ 
nance technique [44] . The Hamiltonian of this problem 
preserves time-reversal symmetry but breaks parity. The 
wave functions of the different spin components are com¬ 
plex in general. In the anti-ferromagnetic domain only 
miscible density profile of different components are found. 
In this domain the wave functions preserve time-reversal 


symmetry. In the ferromagnetic domain phase-separated 
density profile of different components are found. The 
underlying wave functions could be degenerate and break 
time reversal symmetry in this domain with time-reversal 
operator connecting two degenerate wave functions. A 
class of parity-breaking states are found where the real 
and imaginary parts of wave functions exhibit opposite 
parities. These conclusions were illustrated by a numer¬ 
ical solution of a mean-held model 
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